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ABSTRACT 

Standard ACDM predicts that the major merger rate of galaxy-size dark matter halos rises rapidly with 
increasing redshift. The average number of close companions per galaxy, N^ is often used to infer the galaxy 
merger rate. Recent observational studies suggest that Nc evolves very little with redshift, in apparent conflict 
with theoretical expectations for the merger rate of dark matter halos. We use a "hybrid" N-body simulation 
plus analytic substructure model to build theoretical galaxy catalogs and measure Nc in the same way that it is 
measured in observed galaxy samples. When we identify dark matter subhalos with galaxies, we show that the 
observed lack of close pair count evolution arises naturally. This unexpected result is caused by the fact that 
there are multiple subhalos (galaxies) per host dark matter halo, and observed pairs tend to reside in massive 
halos that contain several galaxies. There are fewer massive host halos at early times and this dearth of galaxy 
groups at high redshift compensates for the fact that the merger rate per host halo is higher in the past. We 
compare our results to z '--^ 1 DEEP2 data and to z ~ data that we have compiled from the SSRS2 and the 
UZC redshift surveys. The observed close companion counts match our simulation predictions well, provided 
that we assume a monotonic mapping between galaxy luminosity and the maximum circular velocity of each 
subhalo in our catalogs at the time when it was first accreted onto its host halo. This strongly suggests that 
satellite galaxies are significantly more resilient to mass loss than are dissipationless dark matter subhalos. 
Finally, we argue that while Nc does not provide a direct measure of the merger rate of host dark matter halos, 
it offers a powerful means to constrain both the Halo Occupation Distribution and the spatial distribution of 
galaxies within dark matter halos. Interpreted in this way, close pair counts provide a useful test of galaxy 
formation processes on ^10-100 kpc scales. 

Subject headings: cosmology: theory, large-scale structure of universe — galaxies: formation, evolution, high- 
redshift, interactions, statistics 



1. INTRODUCTION 

Close pairs of galaxies have provided an important empir- 
ical platform for evaluating theories of galaxy and structure 
formation since the early studies of Holmberg ( 1937). Within 
the modern, hierarchical picture of cosmological structure for- 
matio n, galaxy mergers are expected to be relatively common 
(e.g. Blumenthal et al.'1984VLacev & Cole'1993VKola tt eTaP 
[l?99; Gottlober et al. 2001; Mailer etal. 2005) and the enu- 
meration of close galaxy pairs or morphologically disturbed 
systems has long been used to probe the galaxy merger rate 
jZecf&Koo 1989; Burkev et al. 1994; Carlberg et al. 1993 
Woods et a l.iil995[iY ee & Eflingson 1995; Patton et al. iil997t 
Neuschaefe retalJ|19 97: Carlberg et al. 2000; Le Fevr eetaLj 
2000; Patton etal. 2002; Conselice et al. 2003; Bundv et al. 
2004: .Masiedi et aL.2005: .Bell et aL.2006: .Lotz et al.,.2006.) . 



Comparisons of this kind are clearly important, as observed 
correlations between galaxy characteristics and their environ- 
ments suggest that interactions pla y an essential role in set- 
ting many galaxy properties (e.g. Toomre & Toomre 1972; 
[Larson & Tinslevl 119781; iDressW 198 0: Postman & Geller 
Il984; Barton et al. 2000; Barton Giflesp ie et alJ2003h . 

While traditional merger rate estimates have provided em- 
pirical tests of general aspects of the formation of galaxies, 
they have been less useful in testing specific models of galaxy 
formation or in constraining the background cosmological 
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model. One of the primary difficulties lies in connecting the- 
oretical predictions with observational data. For example, the 
standard approach uses the observed close-pair count (or sim- 
ilarly, an observed morphologically-disturbed galaxy count) 
to derive a merger rate through approximate merger timescale 
arguments. This inferred merger rate is then compared to the- 
oretical expectations fox field dark matter halo merger rates, 
which themselves are quite sensitive to the precis e masses and 
mass ratios of halos under consideration ( see ' Mailer et alJ 
2005; Bell et al. 2006, for a discussion of predicted galaxy 
merger rates). In this paper we exploit the information con- 
tained in close-companion counts by adopting a qualitatively 
different approach. We use a "hybrid" analytic plus numeri- 
cal A^-body model to predict close-companion counts directly 
and to compare these predictions with observed companion 
counts. Encouragingly, the results from this method repro- 
duce observed trends with redshift and number density in 
galaxy companion counts. More generally, our results suggest 
that while close companion counts are only indirectly related 
to dark matter halo merger rates, they may be used directly 
to constrain the number of galaxy pairs per host dark mat- 
ter halo. In this way, companion counts provide an important 
general constraint on the occupation of halos by galaxies and, 
thereby, on galaxy formation models. 

In the context of the standard, hierarchical paradigm 
(ACDM ), the merger rates of distinct dark matter halos 
can be predicted robustly. In particular, several studies have 
demonstrated that the predicted rate of major mergers of 
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galaxy-sized cold dark matter halos increases with redshift as 
(1 + z)"', where the expone nt Ues in the ra nge 2.5 ^ m^ 3.5 
(e.g. 'Governato et al.lll99 9': 'Gottlober et al 2001*). Naively, 
one might expect that the fraction of galaxies with close com- 
panions (or the fraction that is morphologically identified to 
be interacting) should increase with redshift according to a 
similar scaling. As we discuss below, the connection between 
merger rates and the redshift dependence of close companion 
counts is not so straightforward. 

Observational analyses using very different techniques to 
measure the redshift evolution in the fraction of galaxies in 
pairs yield a ve ry broad range of evolutionary exponents , 
w ~ - 4 (e.g., 'Carlberg et al.''200^, 'Le Fevre et al.''200^, 
iPatton et al. 2002; Conselice et al. 2003; Bundy et al. 2004). 
With the exception of the Second Southern Sky Redshift Sur- 
vey (SSRS2) at low redshift ( Patton et al-2000.) . these mea- 
surements use incomplete redshift surveys that are often defi- 
cient in close pairs because of mechanical spectrograph con- 
straints. Thus, apparent discrepancies in these studies may 
result from differing definitions of the pair fraction, cosmic 
variance, survey size and selection, and survey completeness. 

One of the most recent explorations of the redshift evolu- 
tion of close-companion counts was performed by the Deep 
Extragalactic Evolutionary Probe 2 (DEEP2) team JLin et alJ 
|2i)04), who reported surprisingly weak evolution in the pair 
fraction of galaxies out to z ~ 1.1. In particular, they reported 
that the average number of companions per galaxy, A'^ , grew 
with redshift as Nc oc (1 +z)'", where m = 0.51 ±0.28. At face 
value, this weak evolution seems to be in drastic conflict with 
the strong, m ^ 3 evolution predicted for dark matter halos. 
Other recent analyses using different data sets an d different 
techriiques yield simil arly "weak" evolution (e.g. IB ell et all 
l2()()6HLotzet alJ200(a) . 

The resolution of this apparent conflict involves the differ- 
ence between distinct halos and subhalos. The "predicted" 
merger rate of m ^^ 3 applies only to distinct dark matter ha- 
los. The prediction for galaxies, on the other hand, is more 
complicated. Distinct halos are predicted to contain ^ 10% 
of their mass in se lf-bound subs tructures known as dark mat- 
ter subhalos (e.g. Klypinetal. 1999). These subhalos are 
the natural sites for satellite galaxies around luminous cen- 
tral objects. Subhalos can orbit within their host halo for 
a time that depends on the satellite-to-host mass ratio and 
on specific orbital parameters, and can appear as close com- 
panions in projection without necessarily signaling an im- 
minent galaxy-galaxy merger Beyond this, precise predic- 
tions for the merger rates of galaxies are sensitive to the 
poorly -understood process of galaxy formation inside dark 
matter halos. In work that is based on dissipationless cos- 
mological simulations, galaxy formation unknowns can be 
absorbed into various prescriptions for assigning observable 
galaxies to dark ma tter halos and subhalos (for related dis- 
cussions see, e.g., 'Klvpin et al."1999^, 'Bullock et al."2000^ 
Diemand et al. 2004; Gao et al. 2004; Nagai & Kravtsov 
2p05; Conrovetal. 2006; Faltenbach er & Diemandl 12006 : 
Wang et al. 2006; Weinberg et al. 2006). 

In this work we use a large cosmological N-body simula- 
tion to calculate the background dark matter halo distribu- 
tion and an analytic formalism to predict the subhalo pop- 
ulations within each of these host halos. There are several 
noteworthy advantages to adopting this approach. First, we 
rely on straightforward mappings between galaxies and dark 
matter halos and subhalos. Second, we measure the average 



number of close companions, Nc, in our simulations in ex- 
actly the same way in which they are observed in the real 
universe. We then use the observed companion counts di- 
rectly, rather than inferred merger rates, to discriminate be- 
tween simple and physically-reasonable scenarios for associ- 
ating dark matter halos and subhalos with luminous galax- 
ies. Our approach also allows us to overcome several techni- 
cal issues. Using analytic models with no inherent resolution 
limits to model substructure allows us to overcome possible 
issues of n umerical o vermerging in the very dense environ- 
ments (e.g.. lKlvpin et al. 1999) where close galaxy pairs are 
often found. Furthermore, our analytic subhalo model allows 
us to assess the variance in close-companion counts that can 
be associated with the variation in substructure populations 
with fixed large-scale structure. Our methods also allow us 
to test for the importance of cosmic variance and chance pro- 
jections. We find these effects to be small for typical close 
companion criteria at redshifts less than z ~ 2 in our preferred 
model. 

A powerful way to quantify and parameterize the way 
in which close pairs constrain the relationship between ha- 
los and galaxies is through the halo occupation distribution 
(HOD) of galaxies within dark matter halos. The HOD is 
the probability that a distinct, host halo of mass M con- 
tains A^ observed galaxies, P{N\M), and is usually param- 
eterized in a simple, yet physically-motivated way (e.g., 
Berlind & Weinberg 2001, Kravtsov et al. 2004a and refer- 
ences therein). Coupling the HOD with a prescription for the 
spatial distribution of galaxies within their host halos allows 
for an approximate, analytic calculation of close pair statistics 
(Bullock, Wechsler, & Somerville 2002). Using pair counts 
to constrain the HOD and the distribution of galaxies within 
halos provides a direct means to separate secure cosmological 
predictions for dark matter halo counts from uncertainties in 
the more complicated physics operating on the scales of in- 
dividual halos and subhalos (e.g., star formation triggers and 
orbital evolution). As we discuss in § 4, the HOD methodol- 
ogy also provides a useful physical platform for interpreting 
our specific simulation results and for extending them to pro- 
vide a general constraint on galaxy formation models. 

The outline of this paper is as follows. In § |2l we out- 
line our methods, discuss our simulations, and investigate the 
importance of interlopers in common estimators of the close 
pair fraction. We present our predictions for the companion 
fraction, Nc, in §|3|and investigate its dependence on galaxy 
number density and redshift. We compare our general model 
predictions to the observed evolution in A^^. with redshift in 
§ 3.1 and with our own analysis of pair counts from the UZC 
redshift survey and SSRS2 at z = in § 3.2. We discuss the 
importance of cosmic variance in pair analyses in §3.3 and 
provide general predictions for future surveys in § 3.4. In §|4] 
we present a detailed discussion of our results in the context 
of the galaxy HOD, and explain the predicted evolution in pair 
counts as a consequence of the convolution of a halo merger 
rate that increases with redshift and a halo mass function that 
decreases with redshift. We conclude in §|5]with a summary 
of our primary results and a discussion of future directions. 

2. METHODS 

The following subsections detail our theoretical methods. 
Briefly, we investigate pair count statistics using an A^-body 
simulation to account for large-scale structure and the host 
dark matter halo population as described in S [ 2.11 We use 
an analytic substructure model JZentner et alJl2005h to iden- 
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Fig. 1 . — The cumulative mass function of host halos as derived from our 
120 /i^'Mpc simulation box at z = 0, 1, and 3. Error bars estimate cosmic 
variance using jackknife errors from the eight octants of the computational 
volume. 



ti fy sa tellite galaxies within these host halos as described in 
§ 12.21 This "hybrid" technique overcomes difficulties associ- 
ated with numerical overmerging on small scales in the large 
cosmological simulation box and has already been demon- 
strated to model accur ately the two-point c lustering statistics 
of halos and subhalos (IZentner et alJl200 y). We use a simple 



method to embed galaxies in our simulation volume (§ 12. 3> 
and "observe" these mock galaxy catalogs in a manner sim- 
ilar to those used in observational studies. We then compute 
pair fraction statistics that mimic those applied to observa- 
tional samples ( I2.4> . 

2.1. Numerical Simulation 

Our simulation was performed using the Adap tive Refine- 
ment Tree (ART) A^-body code JKravtsovet alJ lTQQ?^ for a 
flat Universe with fl,,, = 1 - fi/^ = 0.3, h = 0.7, and cts = 0.9. 
The simulation followed the evolution of 512^^ particles in co- 
moving box of 120/i~'Mpc on a s ide a nd has been previously 
discussed inlAll good et alJJ2006h and lWechsler et alJ J2006I) . 
The corresponding particle mass is mp ~ 1.07 x 10'/j~'Mq. 
The root computational grid was comprised of 5 12-' cells and 
was adaptively refined according to the evolving local density 
field to a maximum of 8 levels. This results in peak spatial 
resolution of /ipeak — 1 .8/;~'kpc in comoving units. 

We identify host halos using a variant of the Bou nd Den- 
sity Maxima algorithm (BDM, 'Klvpin et al.' 199?). Each 
halo is associated with a density peak, identified using the 
density field smoothe d with a 24-particle SPH kernel (see 
iKravtsov et alJ2004aL for details). We define a halo virial ra- 
dius /?vir, as the radius of the sphere, centered on the density 
peak, within which the mean density is Avir(z) times the mean 
density of the universe, pave- The virial overdensity Avii(z), 
is given by the spherical top-hat collapse ap proximation and 
we compute it using the fitting function of B rvan & NormanI 
([1998). In the ACDM cosmology that we adopt for our simu- 
lations, Avh(z = 0) ~ 337 and Avii(z) -^ 178 at z> 1. In what 
follows, we use virial mass to characterize the masses of dis- 
tinct host halos (halos whose centers do not lie within the 
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Fig. 2. — The left-hand vertical axis shows the cumulative number density 
of galaxies in our simulation catalog as a function of velocity using V,ww at 
z = (dot-dash, current Vmax) and V,,, at z = (soHd, the value of Vn,ax when 
the halo was accreted) to identify subhalos as galaxies. The light-hand axis 
maps the z = B-band absolute magnitude of galaxies from the SSRS2 lumi- 
nosity function to the appropriate galaxy number density shown on the left- 
hand vertical axis. Using the right-hand scale, the solid and dot-dashed hues 
demonstrate our adopted z = relationship between Vi„ or Vnow and galaxy 
luminosity. The remaining two lines show the Vi„ sample at z = 1 and z = 3. 
A comparison with Figure 1 reveals dramatic differences between the evolu- 
tion of distinct host halo number densities and the galaxy (subhalo) number 
densities shown here. The errorbars shown were generated by summing in 
quadrature the jackknife error and the realization to realization scatter. 



virial radius of a larger system). The host halo catalogs are 
complete for virial masses M^ 1O"/;~'M0. Figure^shows 
the host halo mass functions at z = 0, 1, and 3 resulting from 
this procedure. The error bars are jackknife errors computed 
from the eight octants of the simulation volume and reflect 
uncertainty in the halo abundance due to cosmic variance. 

2.2. Substructure Model 

In order to determine the substructure properties in each 
host dark matter halo, we model halo mass accretion his- 
tories and track their substructure content using an an- 
alytic technique that incorporates simplifying approxima- 
tions and empirical relations determined via numerical sim- 
ulation. We provide a brief overvi ew of the tec h nique . 
The model is described in detail in |Zentner et al.l (12003 
hereafter Z05), and is an up dated and improved ver- 
sion of the model presented in Zentner & Bullock ( 2003b . 
This model shares many elements with similar models of 
Taylor & Babul (2004), Penarrubia & BensoiJ ( I2005I) . and 
ivan den Bosch et al. c2005iK 

In hierarchical CDM-like models, halos accumulate their 
masses through a series of mergers with smaller objects. Thus 
the first step of any analytic calculation is to build an ap- 
proximate account of this merger hierarchy. For each host 
halo of mass M at redshift z in our numerical simulation vol- 
ume, we randomly generate a mass accre tion history using 
the extended Press-Schechter formalism ("B ond et all 1 199 It 
Lacev & Cole 1993) with the particular implementation of 
ISomerv ille & Kolatt ( 1999). The merger tree consists of a list 
of all of the distinct halo mergers that have occurred during 
the process of accumulating the mass of the final host object 
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at the redshift of observation. Every time there is a merger, 
the smaller object becomes a subhalo of the larger object. 

We then track the evolution of subhalos as they evolve in 
the main system in the following way. At each merger event, 
we assign an initial orbital energy and angular momentum to 
the infalling object according to the probability distributions 
culled from cosmological A^-body simulations by Z05. We 
then integrate the orbit of the subhalo in the potential of the 
main halo from the time of accretion to the epoch of observa- 
tion. We model tidal mass loss with a modified tidal approx- 
imation and dynamical friction using a modifi ed form of the 
Chandrasekhar formula (JChandrasekhail l 1 9431) suggested by 
iHashimoto et al.N200 3'). For simplicity, we model the density 
profiles of halos using the Navarro et al. ( 1997, NFW) profile 
with concentrations set by their accretion histories according 
to the algorithm of Wechsler et al. (2002). As subhalos orbit 
within their hosts, they lose mass and their maximum circular 
velocities decrease as the profiles are heated by tidal interac- 
tions. A practical addition to this evolution algorithm is a cri- 
terion for removing subhalos from catalogs once their bound 
masses become so small as to render them very unlikely to 
host a luminous galaxy. This measure prevents using most 
of the computing time to evolve very small subhalos on very 
tightly bound orbits with very small timesteps. For the pur- 
poses of the present work, we track all subhalos until their 
maximum circular velocities drop below Vmax = 80 km s~', at 
which point they are no longer considered and are removed 
from the subhalo catalogs. We refer the interested reader to 
§ 3 of Z05 for the full details of this model. 

The computational demands of this analytic model are such 
that we can repeat this process several times for each host 
halo in the simulation volume in order to determine the ef- 
fect of the variation in subhalo populations among halos of 
fixed mass on our close pair results. For example, one source 
of noise is due to subhalo orbits: subhalos at pericenter are 
likely to be counted as close pairs, whereas subhalos at apoc- 
enter are less likely to be counted as close pairs. In order to 
make this assessment, we repeat this procedure for computing 
subhalo populations for each host halo in the computational 
volume thereby creating distinct subhalo catalogs with iden- 
tical large-scale structure set by the large-scale structure of 
the simulation. We refer to each of these distinct subhalo cat- 
alogs as a "realization" of the model, a term which derives 
from the inherent stochasticity of the model. Additionally 
we may perform rotations of the simulation volume. These 
rotations provide us with different lines of sight through the 
substructure of the simulation. Since the perpendicular sep- 
aration along the line of sight of the subhalos is extremely 
important to the close pair statistics, along with the velocity 
differences along these lines of sight, these rotations provide 
us with additional effective "realizations". Three rotations are 
performed on each simulation box, providing us with more 
effective "realizations". 

This substructure model has proven remarkably successful 
at reproducing subhalo count statistics, radial distributions, 
and two-point clustering statistics measured in full, high- 
resolution A^-body simulations in regimes where the two tech- 
niques are commensurable. The results of the model agree 
with full numerical treatments over more than 3 orders of 
magnitude in host halo mass and as a function of redshift 
(Z05). This success motivates its use to overcome some of 
the difficulties associated with using a purely numerical treat- 
ment. Specifically, unlike N-body simulations, this method 
suffers from no inherent resolution limits and makes the pro- 



cess of selecting subhalo populations based on specific fea- 
tures of their merger histories or detailed orbital evolution 
clean and easy. Note that any correlation between formation 
history or the HOD and the larger scale environment will not 
be included in catalogs created with this hybrid method, but 
this is not expected to be i mportant for m ost of the mass range 
we consider here (Wechsl er et al.l2006l) . 

We explore two simple yet reasonable toy models for asso- 
ciating galaxies with dark matter subhalos. In order to quan- 
tify the size of subhalos, we use the subhalo maximum circu- 
lar velocity, Vmax = max[y/GM(< r)/r]. This choice is mo- 
tivated by several considerations. As long as only systems 
that are well above any resolution limits are considered, this 
quantity is measured more robustly in simulation data and is 
not subject to the same ambiguity as particular mass defini- 
tions because Vmax is typically achieved well within the tidal 
radius of a subhalo. These facts make our enumerations by 
Vmax easier to compare to the work of other researchers using 
different subhalo identification algorithms. In the next section 
we describe our mapping of galaxies onto subhalos, and ex- 
plore models for mapping galaxies onto subhalos that use the 
maximum circular velocities defined at two different epochs. 

2.3. Assigning Galaxies to Halos and Subhalos 

After computing the properties of halos and subhalos in a 
ACDM cosmology, the next step is to map galaxies onto these 
objects. In our comparisons with observational data, we nor- 
malize our model galaxy catalogs to an observed population 
by matching their mean number densities to those of halos and 
subhalos. In each case, we assume that there is a monotonic 
relationship between halo circular velocity, Vmax, and galaxy 
luminosity in the relevant band of interest. The number den- 
sity of galaxies as a function of halo (subhalo) Vmax is shown 
in Figure|2l For our z = comparisons we adopt the SSRS2 
B-band luminosity function in order to match the samples ap- 
propriately. The implied mapping between Mb and Vmax for 
the SSRS2 luminosity f unction is d emonstrated in the right y- 
axis of Figure|2l For the' Lin et alJ 112004) DEEP2 comparison 
at z ~ 0.5 - 1 . 1 we match our Vmax functions with the measured 
luminosity function in the Rab band. 

Although a monotonic relationship between Vmax and lu- 
minosity is a simplifying assumption, it serves as a useful, 
physical model for our comparisons, and is motivated in part 
by the well-known TuUy-Fisher relation which shows that the 
observed speed of a spiral galaxy correlates with its luminos- 
ity in all bands (Tullv & Fisher 1977). Scatter in this relation- 
ship could (of course) be included in our model, but we have 
chosen to neglect this for the sake of simplicity. Moreover, a 
number of studies have demonstrated that threshold samples 
selected in this way exhibit excellent agreement with many 
of the clustering statistics of observed galaxy samples (e.g^ 
Klvcin et al. 1999; Kravtsov et al. 2004a; Tasitsiomi et all 
12004: .Azzaro et aL.2005: .Conrov et al-2006. ; Marin, Wech- 
sler & Nichol, in preparation). The decision to normalize our 
catalogs based on the number densities of halos and subhalos 
above thresholds in maximum circular velocities circumvents 
the need to model star formation directly in our simulations 
and allows us to focus on well-understood physical param- 
eters. As we demonstrate below, close pair fractions vary 
strongly with galaxy number density at fixed redshift. 

For subhalos, a monotonic relationship between current 
maximum circular velocity and galaxy luminosity is less well- 
motivated than it is for host halos (or halos in the "field"). 
When a halo merges into a larger host halo and becomes a sub- 



Close Galaxy Counts and Hierarchical Structure Formation 



0.3 




0.3 



_ Comoving Coordinates 



3 0.2 

6 



0.1 



1 ' ' ' ' I ' ' ' ' I ' 
Physical Coordinates 




kpc) z 
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as variation in the subhalo populations. The uncertainties were calculated using the realization to realization scatter and the jackknife en'ors. Left: The z = 
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pairs by these observationally-motivated criteria inhabit the same host dark matter halo for z < 1.5. In the physical coordinates we see a sharp rise in the impurity 
beyond this redshift. 



halo, it loses mass as a result of tidal interactions and its max- 
imum circular velocity decreases. The degree to which lu- 
minosity is lost or surface brightness dechnes is considerably 
more uncertain (see, e.g.jBuHock & Johnston 2005, and Bul- 
lock & Johnston, in preparation), but is probably less than the 
decrease in mass (e.g. Nagai & Kravtsov 2005) Using maxi- 
mum circular velocities rather than total bound masses partly 
accounts for this mismatch between mass and luminosity be- 
cause maximum circular velocity scales on ly very slowly 
with bound mass, Vmax ocM^, with 7 ~ 0.25 (iKravtsov et al.l 
l2004bt iKazantzidis et alJ i2004). In addition, environmental 
effects add a further complication as they may influence the 
luminosity or star formation rate of a galaxy in such a way 
that it is driven off of any V^ax -luminosity relationship that 
may exist in the field. 

We deal with this uncertainty by adopting two simple mod- 
els for associating subhalos with galaxies: first, the maximum 
circular velocity that the subhalo had when it was first ac- 
creted into the host halo. Via, and second, the maximum cir- 
cular velocity that the subhalo has at the current epoch, Vnow 
The choice of Vi,, mimics a case where a galaxy is much more 
resistant to reductions in luminosity than the dark matter sub- 
halo is to mass loss. The choice of Vnow represents a case 
where a galaxy drops out of a luminosity threshold sample 
(or, perhaps more realistically, a surface-brightness threshold 
sample) in direct proportion to the decline in the maximum 
circular velocity of the subhalo. 

As we stated above, subhalo Vmax declines much more 
slowly than total bound mass during episodes of mass loss, 
so both methods assume that the luminosity of the galaxy is 
resistant to large changes resulting from the early stages of in- 
teraction, which is supporte d by observations of p airs at low 
redshifts (IBarton et alJ200!llBarton Gillespie etal.120031) . In 



the case of Vin, the luminosity is set in the field and does not 
change after merging into the host system, while in the case of 
Vnow the luminosity declines slowly as a result of interactions 
with the host system. 

The Vin choice is likely the most natural one, and has 
been shown to provide an excellent match to clustering 
statistics over a range of redshifts for scales larger than 
100 /!"'kpc (Conrovetal. 2006). However, the Vnow model 

S irovides acceptable cl ustering statistics at large separations 
Kravtsov et alJl20()4al) . and provides a very useful reference 
point for exploring the uncertainty in theoretical predictions 
of close pair fractions and for highlighting the type of physi- 
cal processes that close pairs can and cannot probe. 

We emphasize that the evolution of the host halo merger 
rate is not an appropriate comparison for the evolution of close 
pair counts, which probe the galaxy or subhalo merger rate. 
The models we examine represent two different choices for 
how galaxies populate their hosts. The fact that the Vin and 
Vnow models both have the same host halo merger rates, but 
different halo occupation counts that result in very different 
predictions for small-scale clustering, highlights the fact that 
close pairs are not directly connected to halo merger rates ex- 
cept in the context of a specific model for how galaxies popu- 
late these halos. 

Figure |2l shows the cumulative number density of "galax- 
ies" identified in our simulations, n^, as a function of their 
maximum circular velocities. The solid line shows z = 
galaxies using Vin as an identifier. The dash-dot line shows 
the same using Vnow- The Vnow choice lies below Vin because 
galaxies more readily fall out of the observational sample in 
this case due to the mass lost by the subhalo after accretion. 
Note that host halos always have Vnow = Vm. The fact that the 
Vnow and Vin number density lines differ by only r^ 30% re- 
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fleets the fact that most galaxies are field galaxies and this 
choice affects only the satellites of group and cluster systems. 
The right y-axis for the solid and dash-dot lines in Figure |2] 
illustrates how this choice affects our mappings between ve- 
locity and galaxy luminosity at z = 0. The dashed and long- 
dash-dot lines show rig for the Vin choice at z = 1 and z = 3. 
The luminosity function used is the combined SSRS2 lumi- 
nosity function jMarzke et al.lll998h . Compared with Figure 
^we see that the galaxy velocity function evolves much more 
weakly with redshift than does the host halo mass function. 

2.4. Defining Close Galaxy Pairs and Diagnosing the Effect 
of Interlopers 

Various definitions have been used fo r iden tifying close 
pairs of galaxies (e g-, Kennicutt & KeeB l 1 9841: IZenf & Kqg 
1^89', 'Burkev et al.'1994; Carlbera et al]ll994t_ IWoods et al 
1^95; Yee & Eflinason 1995; Barton et alJl2000t1P atton et al 
2p00; Carlbera et al. 2000; Bundv et al. 2004; Le Fevre et al 
2000; Patton et al. 2002. 1997; Neuschaefer et al. 1997). In 
this work we follow the definition used in Patton et al. ( 2003) 
and lLin et alJ J2004h . and use the average number of compan- 
ions per galaxy to enumerate close galaxy pairs: 
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A^. 
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(1) 



In Eq. ([3, ng is the number density of galaxies in the sample 
and Hp is the number density of individual pairs. 

We define pairs as galaxies that fall within a well-defined 
line-of-site velocity separation | Av| and that have a projected, 
center- to-center separation on the sky, rp, that lie between an 
inner and outer value r, < rp < ro. The inner radius is chosen 
to prevent morphological confusion. The outer radius may be 
adjusted to reduce contamination of the sample by interlopers. 
We use fiducial values of r,- = 10/i~'kpc and rg = 50/i~'kpc in 
physical units and we adopt a fiducial relative line -of-sight ve- 
locity difference of | Av| = 500 km s"', following Patto n et al.l 
(12002) and Lin et al. (2004). In what follows, we refer to the 
chosen close-pair volume as a "cylinder". We emphasize that 
in our standard measure we adopt physical radii and a con- 
stant Av cut to define our cylinder at all redshifts ("physical 
cylinder"), but we also explore a case where the radii and line- 
of-sight depth (defined by Av) are scaled with the expansion 
of the universe ("comoving cylinder"). 

Our model catalogs provide a useful tool for studying the 
degree of contamination by interlopers in this commonly- 
adopted measure. If we define interlopers as projected galaxy 
pairs that do not lie within the same host halo, we can define 
a measure of close pair count impurity as 



Impurity = — , 

np 



(2) 



where n/ is the number density of interloping, or false pairs 
in the sample and «p is the total number density of observed 
pairs. 

The left panel of Figure |3] shows how the impurity varies 
as a function of the choice of outer pair radius ro for both V,„ 
(solid) and ¥„„». (dashed) samples at z = 0. For each choice, 
we fix r, = 10/;~'kpc and define the galaxy population using 
a galaxy number density ng = O.Ol/i Mpc . The Vnow case 
has a higher overall impurity ('^ 20% compared to r^ 5%), 
reflecting the fact that the amplitude of Nc is generally lower 
for V„ow (see §3 ) and is therefore more affected by chance 
projections. Both the Vin and V„ow samples show only a mild 
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Fig. 4. — Evolution in the companion fraction. Solid sq uares are 
data p oints from the DEEP2 Survey taken from fields 1 and 4, ILin et al] 
2004). Close companions are defined as having physical separations between 
Tp = 1 — 5 /i" ' kpc with a line of sight velocity dif ference of < 500 k m s" ' . 
Solid triangles are data from CNOC2 and SSRS2, IPatton et alj2002 ). The 
empty points with best-fit lines are taken from our simulations. At each red- 
shift we select model galaxies by approximately matching the number densi- 
ties of galaxies used in the DEEP2, SSRS2, and CNOC2 points to the near- 
est simulation redshift. Empty triangles use Vnow velocity function to match 
number densities while empty squares use Vin. Lines show the best-fit to 
Nc (X (1 + z)"' with m = 0.42 ± 0. 17 (V„ow) and m = 0.99 ± 0. 14 (Vin). In both 
cases the predicted evolution is quite weak and broadly consistent with the 
data. 



increase (^ 3%) in the interloper fraction as the outer cylinder 
radius is varied from /•„ = 30 /!"'kpc to 100 /i"'kpc. 

The right panel of Figure |3l shows the impurity as function 
of redshift for the Vin-selected sample. The dashed line corre- 
sponds to a physical cylinder and the solid line corresponds to 
a comoving cylinder at the same fixed comoving number den- 
sity ng = O.Olh^ Mpc"^. Both lines have rg = 50/;~'kpc. The 
impurity in the comoving cylinder remains relatively small 
with z, while for the physical cylinder, the fraction of inter- 
lopers remains small until z '^ 2, after which it rises sharply. 
This reflects the fact that halo virial radii scale as ~ (1 +z)~' 
and become a smaller fraction of a fixed physical r,, at high 
redshift. While one might interpret this result as favoring the 
use of a comoving cylinder for pair identification, there are 
practical difficulties in resolving pairs at small comoving sep- 
arations at high redshift. For example, while physical sepa- 
rations between rp = 10-50 h~^kpc correspond to resolvable 
angular separations at z = 2, 6p ~ 1.7-8.5 arc seconds, the 
corresponding comoving distances would be much more dif- 
ficult to resolve, 9p ~ 0.57-2.8 arc seconds. Ground based 
surveys require separations > 2 arc seconds (Bell et al. 2006). 
The main conclusion to draw from this figure is that for most 
practical pair measures, interlopers will be of increasing im- 
portance at high redshift. 



3. RESULTS 



3.1. 



Evolution in Companion Counts 

We first present results on the overall evolution of close 
companion counts with redshift. Figure ^ shows our model 
results for the companion fraction of galaxies, Nc, as a func- 
tion of redshift for our Vin selection (open squares, with line 
fits) and our Vnow selection criterion (open triangles, with line 
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Fig. 5. — Close companion fraction as a function of number density and redshift for a fixed p/!>'.vica/ separation between r, = lO/i 'kpc and r^ = 50/i 'kpc in 
projection and a Hne-of-sight velocity difference |Av| < 500 km s"'. Small points with linear fits show the the average number of close companions Nc, as a 
function of galaxy number density ng, for our model with galaxy selection based on V„ow (left panel) and Vi„ (right panel). The five lines show power-law fits of 
the companion fraction at z = 0, 0.5. 1, 2, and 3 (from bottom to top) to the foTvciNc oc n" The best-fit slopes vary from a = 0.89 ±0.09 atz = toa = 0.67 ±0.07 
at z = 3 for galaxy samples selected on Vnow and from a = 0.39 ± 0.07 at z = to a = 0.61 ± 0.04 at ^ = 3 for samples selected on thresholds in Vi,,. Large, solid 
points represent observational data points shown in Figure|3 These points are taken from the three nearest redshift bins to z ~ 0, 0.5, and 1.0, and plotted at the 
number densities of the data samples. An'ows along the bottom edge of the figure indicate the redshift corresponding to each of these data points. The decreasing 
number densities with increasing redshift in the observational samples results in a nearly constant value of the companion fraction with redshift. The data seem 
marginally to favor the Vi„ model over Vnow. 



fits) compared with some recent observational measurements. 
Solid squares are data from the Lin et al. (2004) DEEP2 anal- 
ysis (fields 1 and 4). Solid triangles at lower redshift re- 
flect data from the Second Canadian Network for Observa- 
tional Cosmology (CNOC2) survey and data from SSRS2 
JPatton et alJ2002h . 

As was found in the observational sample of iLin et al.l 
(|^04), the companion count exhibits very weak evolution 
when plotted in this manner It is also evident that both the Vin 
and Vnow models show similarly weak evolution and seem to 
be broadly consistent with the observational data. The formal 
fit by Lin et al. (2004) from this data gives close pair evolu- 
tion as Nc oc (l+z)"', with m = 0.51 ±0.28. A similar fit to our 
model results, from z = to z = 1, shown in Figure 0] yields 
m = 0.42 ±0.17 for selection on y„ow, and m = 0.99 ±0.14 
for sample selection based on ¥[„■ Note that the Vin points are 
systematically higher than the Vnow points. This reflects the 
fact that galaxies are more easily driven out of an observable 
sample by tidal mass loss in the Vnow-selected sample (see § 4 
for a more detailed discussion). 

While Figure0]seems relatively simple upon first examina- 
tion, the broad agreement between simulation results and the 
data is attained because we have been careful to normalize the 
data and theory in a consistent manner Each of the DEEP2 
points correspond to different underlying galaxy number den- 
sities. This was done in order to account for gross luminosity 
evolution with redshift, z, assuming a simple model where 
M{z) = M(Q)-z- In order to make a fair comparison with their 
results, we have forced the mean number density of the model 
galaxies to match the corresponding observed number density 
of galaxies by varying the maximum circular velocity thresh- 



old (using either Vin or Vnow)- This technique corresponds 
to re-normalizing the relationship between Vin (or Vnow) and 
galaxy luminosity at each redshift. 

Using the correct number density for comparison at each 
redshift is of critical importance. In Figure |5] we plot the 
companion fraction as a function of comoving galaxy number 
density ng in our model catalogs (small points with line fits) 
using Vnow (left panel) and Vin (right panel) for five redshift 
steps between z = 0.0 and z= 3.0. For both models, at fixed 
z, the companion fraction increases with the overall galaxy 
number density, while at fixed comoving number density the 
companion fraction increases with redshift. The number den- 
sities and close pair counts for selected data points in Figure^] 
are represented as large, solid squares in Figure|5] An approx- 
imate redshift for each of the data points is indicated by an 
arrow along the bottom of the plot. Comparing the lines with 
the solid squares in Figure |5] illustrates that the Ndz) points 
in Figure |4] show little evidence for evolution simply because 
they correspond to lower number densities at higher redshifts. 
An analytic explanation for how A^^ scales with number den- 
sity is given in the Appendix, where we also present z = 
correlation functions for the V^ax models at different values 



of n^. 



3.2. z = Companion Counts 



The trend with A^^ and ng seen in the simulations leads us 
to look for such a trend in some observational samples. To 
compare with z = galaxies, we must balance the need for 
a large comparison redshift survey that fairly samples large- 
scale structure with the requirement for completeness with 
close galaxy pairs. Our approach is to use small but nearly 
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Fig. 6. — The variation of W^ at ^ = based on the number density, «g. 
of galaxies in the sample. We vaiy galaxy number density by making differ- 
ent threshold cuts in V,,, and Vno,,,. Close companions are defined as having 
physical separations between r,, = 10—50 /i^'kpc with a line of sight velocity 
difference of < 500 km s"' . Notice that pair counts based on galaxies selected 
according to V,„ are comparable to data from both the UZC and SSRS2 red- 
shift surveys. In addition, note the radically different slope for the Vnow line 
as well as its under-prediction of Nc . The points along the V,„ line repre- 
sent the average over three rotations of the four realizations at cuts ranging 
from y,„ = 1 10 - 350 in steps of 20 km s" . For Vnow the points range from 
Vnow = 1 10-230 in steps of 20 km s" . The observational data points clearly 
favor the Vm model over the Vnow model. 



complete redshift surveys and to examine the effects of cos- 
mic variance a posteriori. 

In Figure |6l we show the close companion fraction Nc, 
derived from galaxies taken from the UZC (solid squares 
for the North field and solid circles for the South) and the 
SSRS2 and CNOC2 (solid triangles) redshift surveys as a 
function of galaxy number density using the measured lu- 
minosity functions for the UZC-CfA2 North, UZC-CfA2 
South, and SSRS2 s urveys JM arzke et al. 1994; Falco et al. 
[1999; Fasano 1984; Patton et al. 2000). The CfA2 North 
field originally covered a range of declination from 8.5° < 
(5 < 44.5°. For it's Northern fie ld CfA2 has a r i ght as - 
cension range of 8'' < a < 17\ iGeller & Huchral Jl989t) : 
Euchraetal. (1990, 1995). The original CfA2 South field 
covers a range of declination from -2.5° < (5 < 48°. The 
CfA2 survey's southern field has a right ascension range of 
20'' < a < 4\ lO iovanelli & Havnes ( 1985); GiovanelH et d] 
([1986^; 'Havnes et al. (1988); GiovanelH & Havnesi (1198^ 
IWegner et al. ( 1993); GiovanelH & Havnes ( 1993). The CfA2 
North field has 6500 galaxies while the South field has 4283 
galaxies originally cataloged. These are volume-limited sam- 
ples from 2300 < cz < HolO'^^""""' km s"'. Here the ex- 
ponent G(MB.iim) = (15.5-MB.iim-25)/5 is a function of a 
variable limiting magnitude and Mb = 15.5 is the limiting ap- 
parent magnitude of the surveys. We derive number densi- 
ties for each point by integrating the luminosity function from 
MB.iim = -18,-19, and -20, respectively to a minimum mag- 
nitude of Mb = -30. Note that within the same survey the data 
points are not independent. Differences among the surveys 
arise, in part, because of differences in large scale structure, 
see §3.3. 

The results of our model pair counts are also shown in Fig- 
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Fig. 7. — Fraction of unique pair-identified galaxies with agiven lookback 
accretion time. We include each galaxy identified as part of a unique z = 
pair and use = 0.01 tf Mpc"^^ with Vin (solid) and Vnow (dashed). If a pair 
of galaxies sits in different host halos, then each galaxy is assigned an (un- 
physical) negative accretion time and placed in the far left bin (interlopers). 
If a pair of galaxies sits in the same host halo, then each galaxy is assigned 
the time it was accreted into the host. If one of the same-halo paired galaxies 
happens to be a central object, we set the lookback accretion time to be the 
age of the universe (far right bin). 



ure|6lfor both selection on Vin (open squares) and selection on 
Vnow (open triangles). The data clearly favor the higher pair 
counts that follow from selecting galaxies based on the ini- 
tial maximum circular velocities of their subhalos when they 
entered their hosts, Vin- The V„ow selection under-predicts the 
data by a factor of ~ 3 - 5, depending on the number density. 
We emphasize that both initial circular velocity and final cir- 
cular velocity selections are made from the same underlying 
population of host dark matter halos, which thus have iden- 
tical merger histories and instantaneous merger rates. The 
differences in pair counts reflect only differences in the evo- 
lution of satellite galaxies within their host dark matter halos. 
While the host dark matter halo merger rates are identical 
between the V-m and Vnow models, the accretion times of the 
galaxy-identified subhalos in the two cases are significantly 
different. Figure shows the distribution of lookback accre- 
tion times for each object identified in a unique galaxy pair at 
z = for «(, = O.Ol/j'-^Mpc-' using Vin (solid) and Vnow (dashed). 
In this figure every unique galaxy pair is assigned two accre- 
tion times, one for each galaxy involved in the pair. For exam- 
ple, if there are A^ close companions then we include A^(A^ - 1 ) 
accretion times for those galaxies. If a pair of galaxies is con- 
tained in the same host halo, we assign each galaxy a time 
based on its time of accretion into the host halo. If one of 
the same-halo objects is a central galaxy, we let its "accretion 
time" be the age of the universe (far right bin). If the two 
galaxies identified as a unique pair are part of different halos 
(interlopers) we assign each galaxy an unphysical, negative 
accretion time (the far left bin). ^ 

^ Note that if a galaxy is pait of more than one unique pair, its accretion 
time can be included in this diagram more than once. For example, consider 
an object with two close companions, one of which is an interloper in a dif- 
ferent host halo and the other is contained within the same host. We record 
a negative accretion time for this galaxy when we are counting the interloper 
pair and we record a non-negative accretion time when counting it as part of 
the same-halo pair. 
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The far right bin in Figure Q shows that ^ 28% of galax- 
ies in close companions correspond to central objects. How- 
ever, note that if we were enumerating pairs rather then the 
galaxies in the pairs, most pairs on these scales, in our case 
~ 56%, would consist of a satellite and a central object, as 
one might expect from the familiar "one halo term" in the 
halo model. Similarly, ~ 5% (^ 20) of paired galaxies are 
"interlopers" sitting in different host halos (far left bin) for 
the Vin (Vnow) choice. As seen in the central histogram, most 
unique pair-identified galaxies are subhalos. The Vnow-chosen 
satellite pairs tend to be accreted fairly recently, within the last 
~ 0.5-2 /i"' Gyr. On the other hand, the ¥[„ galaxies can sur- 
vive longer within the host potential without dropping out of 
the observable sample, and have a fairly broad range of accre- 
tion times ('^ 1 -4 /i"' Gyr). More generally, we see that the 
inferred accretion time distribution for observed galaxy pairs 
will depend sensitively on galaxy formation assumptions. 

3.3. Cosmic Variance 

Cosmic variance is a major concern for the result presented 
in Figure|6l The UZC and SSRS2 surveys both contain struc- 
tures that are comparable in size to the surveys as a whole, and 
are therefore not "fair" samples of the large-scale structure of 
the universe. For example, the data clearly reflect the fact 
that the SSRS2 survey is underdense relative to UZC. Here, 
we investigate whether the high observed pair fractions could 
simply reflect an overdensity on a scale comparable to that of 
the surveys. 

To examine the uncertainty in pair counts due to cosmic 
variance in more detail we performed mock observations of 
subsets of our full simulation volume to determine the varia- 
tion in the pair counts in volumes comparable in size to the 
survey volumes. Figure |S] shows the variation in Nc over 
eight cubic regions of the simulation for each realization, each 
60/i"'Mpc on a side for a volume of 2.16 x \Q^h~^ Mpc"*. This 
sub-volume is within 12% of the the volume-limited UZC 
survey to Mb iim = -19. It is roughly half of the volume for 
the MBjim = -20 points for UZC-CfA2 North and UZC-CfA2 
South; and more than twice the volume for the MB.iim = -19 
SSRS2 point. For each octant, we plot twelve values of A^^: 
each of the three projections through the box along with four 
statistical subhalo realizations for each projection. The four 
symbol types correspond to each of the four realizations and 
we have introduced a slight horizontal shift to reflect the three 
sets of rotational projections. We select model galaxies by 
fixing the number density to Ug = 0.01/z^ Mpc~^ within the full 
box. This corresponds to a limit of Vin ^162 km s~'. The 
reason for this approach is that it mimics what is done in the 
observational sample, and uses the same effective luminos- 
ity function for the whole volume of the box. If we instead 
choose the same number density within each sub-volume of 
the box, the overall octant-to-octant scatter is reduced by a 
factor of ^ 2. The solid line represents the mean over all 96 
points, while the dot-dashed line shows the RMS variation: 
Nc = 0.063 ±0. Oil. Using only one projection yields a nearly 
identical result. Error bars on each point are Poisson errors 
on the number of pairs and reflect the expected variation in 
pair counts in the absence of cosmic variance. Notice that the 
scatter is somewhat larger than what would be expected from 
shot noise. 

^ Moreover, if we identify 3-d pairs within 50 kpc at the same number 
density (rather than in projection) the fraction of pairs involving a central 
galaxy grows to ~ 76%. 
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Fig. 8. — Testing the importance of cosmic variance in the observational 
data sets. We show Nc measured in each octant (1-8 on the horizontal axis) 
of the simulation volume for each of the four Monte-Carlo subhalo realiza- 
tions (four symbol types at each octant number) and for three rotations of 
each simulation box, the points offset from each octant. En'or bai's on sub- 
halo realizations reflect Poisson errors on the number of pairs. Here we select 
model galaxies using a fixed Vi,, ~ 162 km s"', which con'esponds to a num- 
ber density within the full box of «g = 0.0 1/r^ Mpc"'. The mean, over all 
of the companions fractions, is given by the solid line while the dot-dashed 
lines show the RMS variation, 0.063 ± 0.011. This mean value does not sig- 
nificantly differ from those produced using only one of the possible lines of 
sight. 



As we see from Figure|8l at maximum the octant-to-octant 
(large-scale structure) variation is ^ 50%. This is much 
smaller than the ~ 300% variation required to reconcile the 
Vnow model with the data. We also note that the scatter within 
the local surveys is much smaller than the variation required 
to reconcile the Vnow model with the observational data sets. 
Nonetheless, repeating this experiment with a larger survey 
with better completeness will be required to verify this agree- 
ment. 

3.4. General Predictions 

Figure |5] illustrated our expectation that the average com- 
panion number Nc, within a fixed physical cylinder volume 
increases with redshift at fixed mean co-moving number den- 
sity. However, this increase does not reflect the increased 
merger activity associated with hierarchical structure forma- 
tion. In fact, the predicted evolution is driven primarily by the 
choice of a fixed physical cylinder to define pairs rather than 
a selection volume that is fixed in comoving units. Figure |5] 
shows the identical statistic computed from our model galaxy 
populations using the selection on Vin at z = 0, 1 , and 3, but this 
time using the comoving cylinder (with volume oc {l+zf'). 
The virialized radius of host halos at fixed mass decreases 
with redshift roughly as /?vir fx (1 +z). With these new coor- 
dinates we continue to probe the same "fraction" of each halo 
(of fixed mass) over the entire redshift range. First, notice 
that in comoving units the variation in pair counts with red- 
shift is much more mild than it is in physical units. The counts 
in Figure |5] span at most a factor of two while the counts in 
Figure |5] span more than an order of magnitude over an iden- 
tical redshift range. Second, after normalizing out the effect 
of the cosmological expansion, the pair count per galaxy ac- 
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Fig. 9. — The average companion count as a function of (Vi„-identified) 
comoving volume, iig, using a cylinder with fixed comoving volume for com- 
panion identification. We see that once the expansion of the universe is ac- 
counted for the companion count per galaxy is predicted to decrease with 
redshift at fixed comoving number density (compai'e to Figurelsl which uses 
a fiy^A physical cylinder volume in defining Nc)- 



tually decreases with redshift rather than increasing as shown 
in Figure |5] 

The primary reason why the companion fraction does not 
evolve strongly with redshift is that the number density of host 
halos massive enough to host more than one bright galaxy de- 
creases at high redshift. This is shown explicitly in Figure FTOI 
where we plot the distribution of host halo masses for galax- 
ies (solid lines with error bars, left axis) at z = and 3 using 
Vin to fix rig = O.Ol/i^^Mpc"^ . The dashed line in each panel 
(right vertical axis) shows the distribution of host halo masses 
for pairs identified in the same sample of galaxies. At all red- 
shifts, pairs are biased to sit in the most massive host halos. At 
high redshift massive halos become rare, and the pair fraction 
is reduced accordingly. 

Given that halos with multiple bright galaxies are expected 
to become increasingly rare at high redshift, faint companion 
counts may provide a more useful probe. Figure^Jshows the 
evolution in the average number of "Faint" companions per 
"Bright" galaxy as a function of redshift using our standard 
fixed physical separation criteria for defining close pairs. The 
top (dot-dashed) line shows the average number of "Faint" 
close companions with Vin > 120 km s"' (Mb < —17.3 at z = 0) 
around "Bright" galaxies with Vi„ > 200 km s"'(Mb < -18.9 
at z = 0) as a function of z. The four lines below this re- 
strict companions to be "brighter" galaxies with Vin = 140, 
160, 180, and finally 200 km s"'. The last of these four 
(solid line) represents our standard "Bright-Bright" compan- 
ion count, Nc- Finally, the lowest line shows the evolu- 
tion in "Very Bright- Very Bright" companion counts using 
Vin = 280 km s"' (marked "280/280"). As expected, the faint 
companion fraction is systematically higher than the bright 
companion fraction, but it also evolves more strongly with 
redshift. On average, almost every "Bright galaxy" (with 
Vin > 200 km s"') at z > 2 is expected to host a "faint" com- 
panion (with Vin > 120 km s"'). The "200/120" Bright-Faint 
companion count rises by almost a factor of '-^ 10 between 
z = and z = 2, while the increase is only a factor of '^ 3 



for the 200/200 case. Indeed the brightest "Bright-Bright" 
companion count (280/280) actually declines beyond z^ 1 .5. 
These trends reflect both the decline of massive halos at high 
redshift and the competition between accretion and disruption 
of subhalos, which favors accretion at high redshift. 

4. INTERPRETATIONS AND THE HALO MODEL 

As demonstrated above, predicted (physical) close pair 
counts do not evolve rapidly with redshift, even at fixed global 
number density. This result is driven by a competition be- 
tween increasing merger rates and decreasing massive halo 
counts: while the number of galaxy pairs per host halo in- 
creases with z (as the merger rate increases), the number of 
halos massive enough to host a pair decreases with z. In this 
section we will work towards a more precise explanation and 
compare our two models for the connection between galaxy 
light and the subhalo properties Vin and Vnow, by investigating 
their Halo Occupation Distributions (HOD). 

The HOD is the probabihty, P(N\M), that A^ galaxies meet- 
ing some specified selection criteria reside within the virial 
radius of a host dark matter halo of mass M. In what fol- 
lows, we provide a brief introduction to the HOD as used in 
galaxy clustering predictions and use a "toy" HOD model to 
show that galaxy close companion counts provide an impor- 
tant constraint on the HOD and the distribution of galaxies 
within halos (see also Bullock et al. 2002). We then quantify 
our two subhalo models, Vin vs. Vnow, in terms of these inputs 
and present this as a general constraint on the kind of HODs 
required to match the current pair count data. 

4.L Pair counts via the halo model 

The halo model is a framework for calculating the cluster- 
ing statistics of galaxies by assuming that all galaxies lie in 
dark matter halos (e.g . Seljak 2000; Peacock & Smith 200(J 
IScoccimarro et al.l2001i) . This approach relies on the fact that 
the clustering properties and number densities of host dark 
matter halos can be predicted accurately, and uses a param- 
eterized HOD to make the connection between dark matter 
halos and galaxies. 

Consider an HOD which consists of a "central" and 
"satellite" populatio n of galaxies orbiting in the same 
dark matter halo (e.g. Kauffmann et al. 1993;'Kravtsov et alj 
2004a; Tinker etal. 2005; Cooray & Milosavljevic 2005). 
The total HOD is then described by two distributions 
P{Nc\M) and P(Ns\M) with individual mean values NdM) = 
J NcP{Nc\M)dNc and Ns(M) = J NsP{N,\M)dN, respectively. 
Let us assume that we can write the first moment of P(N\M) 
(the average number of galaxies per host halo of mass M) as 



N{M)=Nc{M)+N,{M) = 




M>Mi) 
M <M(), 



(3) 



where Mq is the minimum host mass large enough to host 
an observable galaxy. Mi is the typical host mass which 
hosts one observable satellite galaxy, and p is a power which 
describes the scaling of satellite galaxy number with in- 
creasing host mass. Though simplified, this general form 
with p sa 1 is motivated by both analytic expectations (e.g 
IWechsler et^ljj2(^ ^ Z05) and the expectations o f simula- 
tions re.g.lKravtsov et alJ2004al:lTinker et alJ2()05l). and pro- 
duces large-scale clustering results that are in broad agree- 
ment w ith data selected ov er a range of luminosities and red- 
shifts (Conrov et al. 2006). These results also motivate us to 
assume that the satellite HOD is given by a Poisson distribu- 
tion and that the central piece is represented by a sharp step 
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Fig. 1 1 . — The average number of close companions per bright galaxy 
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the brightest "Bright-Bright" companion count shows a definite drop beyond 
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function (in this approximation) at M > Mq. The left panel of 
FigureElshows a cartoon representation of this HOD. 

With P(N\M) given, the number density of galaxies can be 
written as an integral of N(M) over halo mass, 



Mo 



dM 



(4) 



where, dn/dM is the host dark matter halo mass function. As 
shown in Figure^ the host halo mass function is a rapidly de- 
clining function of mass. This implies that the overall galaxy 
number density is dominated by "field" galaxies in host halos 
of mass Ml < M< Mq, where N{M) ~ 1 . 

On the other hand, for small separations, the number den- 
sity of close pairs, Up, will be dominated by pairs contained 
within single halos (Figure|3lillustrates that, for the V-m model, 
^ 10% of close pairs reside in separate host halos according 
to typical close-pair definitions out to z < 3). Ignoring for 
the moment the precise velocity selection criteria, we can ap- 
proximate the number density of pairs with physical projected 
separations meeting some range ri < r < r2 as 

1 /"OO 7 

«p~«pah=-y dM^{mN-\))^F{ri,2\M), (5) 



where {N(N- 1))^ is the second moment of the HOD and 
F(ri,2\M) is the fraction of galaxy pairs that have projected 
separations between ri and r2 within a dark matter halo of 
mass M. The function F provides an ^ order unity covering 
factor determined by the distribution of galaxies within ha- 
los 4. Notice that since (A^(A^- 1))^ = N,{Mf and NAM) -^ 
for M <M\, the pair number density Up (unlike Ug) explicitly 

■* More precisely, F is determined by convolving the projected density dis- 
tribution of galaxies within halos with itself. 
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Fig. 12. — l^ft: A cartoon representation of an input halo model N{M), where the plateau at N{M) = 1 represents central objects and the power-law piece, 
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with rig is steeper because satellite galaxies tend to be more plentiful in lower-mass host halos. The line labeled "c = 5" has a more diffusely packed satellite 
distribution and as a result produces fewer close pairs with all other parameters held fixed. 



selects against the mass regime dominated by "central galax- 
ies" and favors halo with M > M\. 

This simple and intuitive formalism immediately provides 
insight into the nature of the companion fraction, Nc = 
In^/ug. In the right-hand panel of Figure^] we show the 
predicted value of A^^. vs. Ug as calculated using the (approxi- 
mate) equations|4]and|5]with r\ = 10/i~'kpc and ra = 50/i~'kpc 
for several choices of HOD parameters defined in Equation 
|3l In all cases we assume that galaxies are distributed in 
dark matter halos following NFW profiles. The upper long- 
dashed line and solid line assume Mi/Mq =10 and 25 re- 
spectively, use /? = 1 in the satellite Ns(M) piece, and as- 
sume the N FW concentration-mass relation, c(M), given in 
iBuUock et al._ (.2001.) . The dot-dashed line is the same as the 
soUd line except we fix the NFW concentration to a relatively 
low value, c = 5, for all halos. Finally, the lower dashed line 
assumes Mi /Mq = 100 with /:> = 1 .2 and the standard concen- 
tration relation. 

From Figure^lwe see that as the size of the "plateau" re- 
gion increases (characterized by increasing the ratio Mi /Mq) 
the galaxy number density is increasingly dominated by cen- 
tral galaxies and the pair fraction goes down. Also, even for a 
fixed HOD, Nc decreases as galaxies are packed together more 
diffusely in halos (characterized by lowering c). Finally, as 
the satellite contribution steepens (higher p), Nc increases at 
higher number densities because lower mass host halos begin 
to contain multiple galaxies. 

As illustrated by this toy-model exploration, pair count 
statistics should provide an important constraint on galaxy 
HODs as well as the distribution of galaxies within halos. 
The constraint could prove particularly powerful when com- 
bined with larger-scale (r^ 100 kpc) two-point correlation 
function constraints on the HOD. Correlation function con- 



straints generally suffer a severe degeneracy between P(N\M) 
and F{r\M), but this degeneracy may be broken when coupled 
with close {r"^ 50 kpc) pair fraction constraints. 

4.2. Discussion of Model Results 

The halo model serves as a useful framework for discussing 
and explaining our simulation model results. FigurelT^shows 
the first and second moments of our simulation model HODs. 
The upper panels show N{M) and the lower panel shows the 
number of pairs per halo compared to the Poisson expecta- 
tion, < A^(A^- 1) >'''^ / < A^ >. The left side of the figure 
shows how the HOD changes as we vary the selection cri- 
teria, with Vin > 120km s"'subhalos represented by the solid 
line and Vnow > 120km s"'subhalos shown by the dashed line. 
As must be the case, the central galaxy term (analogous to Mq 
in our toy example above) is unaffected by this choice, while 
the satellite component (points) is enhanced when galaxies 
are selected using Vin. This has the effect of extending the 
"plateau" region in the N{M) function and is qualitatively sim- 
ilar to increasing the M\ /Mq ratio in our toy-model discussed 
above. In the right hand panel we show N{M) for two dif- 
ferent velocity cuts (120 compared to 200 km s"') both using 
Vin as the velocity of relevance. Here we see that both the 
satellite term and the central term are affected in this case (Mq 
becomes larger as more massive satellites are selected). Note 
that as expected from the previous discussion of satellite vs. 
central galaxy statistics, the average number of pairs of galax- 
ies per halo tracks the Poisson expectation for host masses 
where N(M)^ 2 and becomes substantially sub-Poissonian 

for N(M)^ 2. As discussed in reference to Equations |3 and 
|5J this means that while the total number density of galax- 
ies (with, e.g., Vin^ 120km s"') will be dominated by host 
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Fig. 13. — Top: The average numberof galaxies per host halo as a function of host halo mass. Bottom: Expression for the pairs per halo <A'(A'-1) >''- /A'(M) 
(relative to the Poisson expectation) as a function of host halo mass. In the upper panels, points with error bars reflect satellite galaxies and Hnes show the total 
(satellite plus central) galaxy count per halo. The left figure compares the V[^ model (solid) and Vnow model (dashed) using the same velocity cut to define 
the galaxy population (120 km s"'). The right panel compares the V|„ model at two different velocity cuts (120 and 200 km s"'). The pair count per halo is 
well described by a Poisson distribution for the satellites and substantially sub-Poissonian for the central objects. Pair counts thus become important only for 
N(M)^ 2, and will be dominated by massive halos. 



halos where N(M) ~ 1 or with M ~ Mq ~ 1O"/!"'M0, the to- 
tal number density of close galaxy pairs will be dominated 
by much more massive host halos M^ Mi ~ 5 x 10'^/i~'Mq 
(where N(M)> 2). 

In addition to its effect on P(M\N), whether the galaxy 
light is more accurately traced by Vnow or Vi,, will also af- 
fect the relative distribution of satellite galaxies within host 
halos. Recall that V,iow selects "observable" galaxies based 
on a subhalo's current maximum circular velocity and V-m se- 
lects galaxies based on the subhalos' circular velocity when 
it was accreted into the system (these might mimic models 
where galaxies are more or less affected by mass loss or sur- 
face brightness dimming after accretion). As shown in Figure 
1141 the Vin selection results in a satellite galaxy distribution 
that is more centrally concentrated (and thus has more close 
pairs per halo) than does the Vnow selected sample. The fig- 
ure shows cumulative number of satellites relative to the to- 
tal number of satellites as a function of radius, scaled to the 
virial radii of the hosts, from the host halo center using host 
halo masses of Myir = 1O'^~'^^/!~'M0. The solid line shows 
this for the Vin case and the dashed line reflects the Vnow case. 
For reference, the dot-dashed line above shows the cumula- 



tive mass distribution for a typical dark matter halo of this 
size ( Bullock et al. 2001). The Vnow case shows the effects of 
enhanced satellite destruction near the center. This difference 
certainly affects predictions for pair counts, however the main 
difference between Nc counts in our two models is driven by 
the differences in their HODs. 

Finally we turn to the evolution of the companion fraction 
with redshift. As discussed in the previous section, the aver- 
age close companion count per galaxy is expected to evolve 
only weakly with redshift (e.g. Figure |4} and even declines 
at fixed comoving number density when pairs are identified 
within a fixed comoving separation (Figure |9j. This occurs 
even though the merger rate per host dark matter halo in- 
creases at high redshift. This increase in merger rate is indeed 
reflected in the HOD. As shown in Figure ^J the average 
number of galaxies per halo at fixed host mass increases sys- 
tematically at high redshift. However, from Equation [S] the 
number density of pairs depends both on the number of pairs 
per halo (< A^(A^- 1) >'^ Ns(M)^) and the number density of 
host halos (dn/dM). In Figure^] for example, N{M) is a fac- 
tor of ^ 2 higher at z '^ 3 at the host mass where pair counts 
become non-negligible Mhost ~ 10'^/z~'Mq. However, host 
halos of this size are a factor of ^ 5 rarer at z = 3 compared to 
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Fig. 14. — The normalized cumulative number of satellite galaxies 
as a function of radius from the host halo center for host halos of mass 
10l2-i2.25;,-l^g (R„, ^ 204/i-'kpc). The solid and dash-dot lines corre- 
spond to Vjn and Vnow = 120 km s"' , respectively. The cumulative mass pro- 
file for a typical halo of this mass is shown by the dot-dash line. En'or bars 
are RMS variations over the four Monte-Carlo subhalo realizations. 



z = (Figure^. This results in a mild overall decline in the 
number of same halo pairs at high redshift. 

The HOD formalism also helps us understand Figure [TTl 
where we find that the average number of "Faint" galaxies per 
"Bright" galaxy increases more rapidly than does the "Bright- 
Bright" fraction with redshift. The more pronounced evolu- 
tion comes about because "Faint" satellite selection decreases 
the characteristic host mass required for a satellite (Mi j). 
At the same time, we focus on "Bright" central galaxies, 
which keeps the minimum host halo mass for hosting a cen- 
tral galaxy (Mo) roughly fixed. As Mi /Mo decreases, we are 
more likely to see an increased companion count (Figure[T2t. 

5. CONCLUSIONS 

In this paper we have investigated the use of close pair 
counts as a constraint on hierarchical structure formation 
using a model that combines a large (120/i~'Mpc box) 
ACDM N-body simulation (Allgood et al. 2006) with a rigor- 
ous analytic substructure code (Zentneret al. 2005) in order 
to achieve robust halo identification on small scales within a 
cosmologically-relevant volume. We measure close pairs in 
our simulation in exactly the same way they are observed in 
the real universe and use these close pair counts to test simple, 
yet physically well-motivated models for connecting galaxies 
with dark matter subhalos. We explore two models for the 
connection between galaxies and their associated dark matter 
subhalos: one in which satellite galaxy luminosities trace the 
current maximum circular velocities of their subhalos, Vnow, 
and a second in which a satellite galaxy's luminosity is tightly 
correlated with the maximum circular velocity the subhalo 
had when it was first accreted into a host halo, Vin. The latter 
case corresponds to a model where galaxies are much more 
resilient to mass loss than are their dark matter halos. 

A highlighted summary of this work and our conclusions 
are as follows: 

• We showed that ACDM models naturally reproduce 
the observed "weak" evolution in the close companion 
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Fig. 15. — The average number of galaxies per host dark matter halo with 
Vin > 120km s^'shown at z = 0, 1, and 3. The points reflect satellite halos 
while the lines include all galaxies. 



fraction, A^,., ou t to z ^ 1.1 as reported by the DEEP2 
team (Lin et al1 l2004l) . The result is driven by the fact 
that while merger rates (and the number of galaxy pairs 
per halo) increase towards high redshift, the number of 
halos massive enough to host a bright galaxy pair de- 
creases with z- 

• We used the SSRS2 and UZC surveys to derive com- 
panion counts as a function of the underlying galaxy 
population's number density at z = and showed that 
the relatively high companion fraction observed favors 
a model where galaxies are more resistant to mass 
loss than are dark subhalos (~ the Via model). Cos- 
mic variance is unlikely to be qualitatively important in 
this conclusion, although verifying the result requires a 
complete redshift survey with a larger cosmic volume. 

• We argued that the close luminous companion count 
per galaxy {Nc) does not track the distinct dark matter 
halo merger rate. Instead it tracks the luminous galaxy 
merger rate. While a direct connection between the two 
is often assumed, there is a mismatch because multiple 
galaxies may occupy the same host dark matter halo. 
The same arguments apply to morphological identifi- 
cations of merger remnants, which also do not directly 
probe the host dark halo merger rate. 

• We showed that while close pair statistics provide a 
poor direct constraint on halo merger rates, they provide 
an important constraint on the galaxy Halo Occupation 
Distribution (HOD). In this way, pair counts may act as 
a general tool for testing models of galaxy formation on 
< 100 kpc scales. 

Our results are not driven by the effects of false pair identi- 
fication. We showed that standard techniques for identifying 
close pairs via close projected separations and line-of-sight 
velocity differences do quite well in identifying galaxy pairs 
which occupy the same dark matter halos (Figure |3}- Typ- 
ically > 90% of pairs identified in this manner occupy the 
same host dark matter halos out to a redshift of ^^ 2 using the 
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Vin model. It should be noted that this fraction does depend on 
the details of the association between subhalos and galaxies. 
The striking differences between our two toy models provides 
evidence of this. Typically, ^ 20% of pairs occupy different 
halos in the Vnovj model. 

In conclusion, current close pair statistics are unable to rule 
out a simple scenario where galaxies are associated in a one- 
to-one way with dark matter subhalos (in particular our Vin 
association) in the current concordance ACDM model. Bet- 
ter statistics, more complete surveys, and the ability to divide 
samples by color or stellar mass indicators will help refine 
this model and explore galaxy formation unknowns on < 100 
kpc scales. Enumerating/a/nf close companion counts around 
bright galaxies (Figure [TTl will also provide a new and pow- 
erful test. 
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APPENDIX 

PAIR FRACTION AS A FUNCTION OF NUMBER DENSITY 

One of the important findings in our analysis is that the observed companion fraction of galaxies should increase with the 
underlying galaxy number density (see, e.g. Figure|5}. We can understand the tendency for Nc to increase with rtg at fixed redshift 
by considering the number of pairs expected in a fixed (physical) cylinder of volume y (e.g., with length £ ~ 2 Av//q' ~ 10/;"' Mpc 
and radius r ~ 50 kpc). Let us start with the idealization of an uncorrelated sample. For an uncorrected galaxy population, with 
mean number density Ug in a large fiducial volume Vq, the total number of galaxies in the sample is clearly 



A? 1 = n Vb. 



(Al) 



The average number of companions per galaxy is A^^ = UgV, so the number of companions in the entire sample is simply n^VVo. 
Therefore, for an unclustered galaxy population, the average number of companions per galaxy should increase linearly with 
galaxy number density, 

Nc = ngV oc rig. (A2) 

Extending this argument to the case of a clustered galaxy population is straightforward: 

N, = ngV[l+0, (A3) 

where ^ = V"' Jy ^ir)dV is the average of the correlation function over the volume of the cylinder. The total number of compan- 
ions in the sample is then is n^[l + ^]Wb and the companion fraction is 



Nc = ngV[l+0^ngV^, 



(A4) 
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where we assume £,{r) :^> 1 at all relevant scales. There are both theoretical and dkect observational reasons to believe that ^ 
should increase with decreasing number density. If we take a simple power-law description ^ oc n"" then we have 



Nc oc n; 



(A5) 



Therefore, quite generally, as long as a < 1, the companion fraction should be an increasing function of galaxy number density. 
The z = correlation function for our two simulation models at two characteristic galaxy number densities is shown in the left 
(Vnow) and right {Vi„) panels of Figure IAT61 Note that ^(r) at r^ 100 /!"'kpc is quite flat for the Vno^ case but continues to rise 
for the Vin case. This reflects the enhanced "destruction" included in the Vnow model (galaxies in the centers of halos tend to lose 
mass quickly and their Vmax values drop out of the sample.) If we characterize the clustering vs. number density trend by fitting 
the amplitude of the z = correlation functions at 100/i~'kpc with ^loo oc riZ" we find a ~ 0.15 for Vnow and a ci 0.5 for Vin. From 
the above arguments with ^loo ~ ^we would then expect A^^ oc «"^^ and Nc oc rfi-^ for Vnow and Vin respectively. Indeed, these are 
approximately the slopes measured at z = in Figure|5] 
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